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ABSTRACT 



The acoustic power injected by turbulent convection 
into solar-like oscillations depends on the details of 
the turbulent spectrum. A theoretical formulation 
for the oscillation power is developed which general- 
izes previous ones. The formulation is first calibrated 
on a solar model in such a way as to reproduce the 
solar seismic data. This allows to investigate dif- 
ferent assumptions about the stellar turbulent spec- 
trum. We next explore consequences of the assumed 
turbulent description for some potentially solar-like 
oscillating stars. Large differences are found in the 
oscillation power of a given star when using different 
turbulent spectra as well as in a star to star compar- 
ison. Space seismic observations of such stars will be 
valuable for discriminating between several turbulent 
models. 



1. INTRODUCTION 



The acoustic power injected into the p modes by tur- 
bulent convection has been modeled by several au- 
thors (iBalmforth, 1992fc |Goldreich et ail, 1991 |GuT 
|ireich fc Keeley 1977] ). In |Samadi fc Goupil (2000| J 



a general formulation is proposed which allows to 
investigate more consistently different assumptions 
about the stellar turbulence such as its turbulent en- 
ergy spectrum. 

Providing that accurate measurements of the os- 
cillation amplitudes and damping rates are avail- 
able it is possible to evaluate the power injected 
into the modes and thus- by comparison with the 



observations- to constrain current theories (Samad 
|fc Houdek, 2000| ). 



In the present paper, the formulation, viewed as a di- 
agnostic tool of stellar turbulent spectra, is used to 
compute oscillations power of several solar-like stars. 
We show that the expected low detection thres hold of 



a space seismic experimen t such as COROT ( Baglin 
[fc The Corot Team, 199<t ) for instance will provide 
highly accurate acoustic power spectra and therefore 
important information on the associated turbulent 



spectrum if the solar-like oscillating targets are prop- 
erly chosen. 



POWER INJECTED INTO SOLAR-LIKE 
OSCILLATIONS 



2.1. Stochastic excitation 

The acoustic power injected into t he oscillations is 
defined (e.g. Goldreich et al., 1994 ) in terms of the 
damping rate ij, the mean-square amplitude (A 2 ), 
the mode inertia / and oscillation frequency lo : 



P(io) = (A 2 ) I lo 2 



(1) 



The mean-square amplitude accounts for both the 
excitation by turbulent convection and the damping 
processes. It can be written in a simplified expression 
as 



(A 2 ) oc rr 1 



M 



dm pq w 



dr 



§(oj,m) (2) 



where £ r is the radial displacement eigenfunction, po 
the density, w the vertical rms velocity of the con- 
vective elements and S the turbulent source function 
accounting for both the Reynolds and the entropy 
fluctuations. Detailed expressions for (A 2 ') and S 
are given in |Samadi fc Goupil (2000 ). The source 



function involves the turbulent kinetic energy spec- 
trum E(k) and the turbulent spectrum E s (k) of the 
entr opy fluctuations whic h can be related to E(k) 
(e.g. samadi et al., 2000a). The turbulent spectra in 
§ are integrated over all eddy wavenumbers k and S 
is in turn integrated in Eq.(|2]) over the stellar mass 
M. 



2.2. Models for the solar turbulence 

From the turbulence theory it is expected that E{k) 
follows the Kolmogorov spectrum as E(k) cx k p with 
the slope p = —5/3. Observations of the solar gran- 
ulation allow one to determine the turbulent kinetic 



2 



spectrum E(k) of the Sun. Observations of t he so 
lar granulation by | Espagnct ct al. (1993) and jN 



OSl; 



ct jl. (1993|) . confirm the existence of a turbulent 



cascade with p ~ —5/3. On the other hand, in the 
k range (small k) prior to the turbul ent cascade of 
slo pe p ~ — 5/3, [Espagnct ct al. (1993) determined a 



The power P(u>) is computed for a calibra ted so- 
lar m odel obtained with the CESAM code ( Morel 
1997 ). Convection is describ ed according to the clas - 
sical mixing- length theory ( Bohm - Vitense, 1958 ) 



reg me with a slope close to 0.7 whereas Nesis et al 



(19 )3|) found a slope of p = —5. 



The oscillation properties were obta ined from the 
adiabatic F ILOU pulsation code of [Trail Minh"fc| 
Leon (1995|) . The physical ingredients o f the model 



are detailed in Samadi & Houdek (2000) 



The free parameters are then adjusted to obtain the 
best fit of the frequency dependence and the max- 
imum amplit ude to the solar observations by Lib- 




Figure 1. Kinetic turbulent spectra versus wavenum 
her k. 



The spectrum observed by Espagnet et al. (1993) is 
modeled by the so called "liaised Kolmogo rov Spec- 
trum " (RK S hereafter) as suggested by Musiclak 



ct i 
tra: 
after 
et 



. (1994 ). We also consider two additional spec- 
the "Nesis Kolmogorov Spectrum" (NKS here 



in ag reement with the observations of Nesis 
i d. (1993| ) - and the "Broad Kolmogorov Spec- 
trum" (BKS hereafter). These spectra (Figure lj) 
obey the Kolmogorov law for k > ko where ko is the 
wavenumber at which the turbulent cascade begins. 

The wavenumber kg is unknown but can be related 
to the mixing length A as fco = 27r / (j3A) where j3 is 
a free paramete r introduced for the arb itrariness of 
such definition ( Samadi fc Goupil, 2000| ). Moreover 
the definition of the eddy time correlation, which en- 
ters the description of the turbulent excitation, is 
somewhat arbitrary and is therefore gauged by in- 
troducing an additional free parameter A. We show 
that the oscillation power computed on the Sun is 
very sensitive to the value s of the free parameters A, 
(5 in Samadi et al. (2000a ). 



v [mHz] 

Figure 2. Computed surface velocity v s (Eq$) as- 
suming the different turbulent spectra of Fig. (0). A 
and (3 values result from fittin g the computed to the 
solar seismic observations by Libbrecht (198t ) using 
the observed damping rate n. 

Results of the fitting are shown in Fig. 0. All spec- 
tra of Fig fit well the solar observations at low fre- 
quency iy < 3.5 mHz) while the main differences are 
observed at high frequency. The overall best agree- 
ment is obtained with the NKS. Details of the result- 



ing adjustments are given in Samadi et al. (2000a). 



2.3. 



The solar case: comparison with observations 
and calibration of the free parameters 



The power injected into the solar oscillations is re- 
lated to the rms value v s of the surface velocity as 



vl = g{r s )P/2riI 



(3) 



where r s is the radius at which oscillations are mea- 
sured. Observations of the solar oscillations provide 
v s and 77 such that P(oj) can be evaluated according 
to Eq.(§). 



3. SCANNING THE HR DIAGRAM 



3.1. The models 



We focus on low intermediate mass stars (1 < M < 
2M & ) because the existence of an outer convective 
shell enables stochastic excitation. The thickness 
of the convective shell depends in particular on the 
star luminosity (L) and effective temperature (T e g). 
We thus consider six models well distributed in the 
solar-like oscillation region (Table Q). These models 
are based on the fundamental parameters of several 
bright stars which w ere selected as targets for EVRIS 
( Baglin ct al., 1993] ) . The equilibrium models and the 
associated eigenfunctions of these stars are obtained 
in the same way as for the solar model. 
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Models 


L 


Tog 


M 


v c 


Age 






[L & ] 


[K] 


[M e ] 


[mHz] 


[Gyr] 


A 


a Tri 


12.1 


6350 


1.68 


1.0 


1.79 


B 


T) Boo 


9.0 


6050 


1.44 


1.0 


3.05 


C 


Procyon 


6.6 


6400 


1.46 


1.5 


2.40 


D 


/3 Hyrdi 


3.7 


5740 


1.08 


1.5 


7.33 


E 


(3 Vir 


3.5 


6120 


1.25 


2.3 


4.10 


F 


7r 3 Ori 


2.6 


6420 


1.25 


3.6 


1.76 



Table 1. Stellar Parameters of selected stars. v c is 
the cut-off frequency. For each model we have given 
the corresponding bright star from the EVRIS pre- 
liminary list of targets. 



3.2. Star to star comparison 

The acoustic power is computed for each selected star 
and for each of the three turbulent spectra of Fig.l 
and is plotted for four of them in Fig.[| Of our 5 
stars, model E is the closest to the Sun in the HR 
diagram (Figure |J), hence its structure (assuming 
here the same chemical composition) is nearly the 
same than in the solar case. Accordingly the oscilla- 
tion spectrum of model E is quite similar to the solar 
one (except small differences at high frequency) and 
therefore does not show large differences when using 
different turbulent spectra (due to the solar calibra- 
tion mentioned in Sect. 2.2); only small differences 
can be seen at high frequency. Same conclusion can 
be drawn for model D. 

For a hotter star, the oscillation spectra computed 
with different turbulent spectra differ at low frequen- 
cies. As it is illustrated with the error bars and the 
detection threshold in Fig. 2, such differences are sig- 
nificant enough to constrain stellar turbulent spectra 
providing that accurate measurements of P(to) arc 
available. 

Differences in oscillation power between model E 
(equivalently the Sun) and the other stars is related 
to the more extended outer convective zone (CZ here- 
after) of model E. As the thickness of its CZ is 
larger, excitation of the low frequency oscillations 
involves all turbulent eddies. Therefore as for the 
Sun the power injected into low frequency oscilla- 
tions is mainly governed by the eigenfunction behav- 
ior (d^ r /dr and I) and not by the behavior of the 
turbulence in the energy injection range (i.e at large 
turbulent scale). 

For the other hotter stars, the relative smaller size of 
their CZ causes the excitation region to be localized 
in a thinner domain where changes in the properties 
of £ r are large. At a given frequency the extension 
of the excitation region changes with the considered 
turbulent spectra. Therefore the shape of P(w) is 
very sensitive to both the properties of £ r - thus to 
oj - and the behavior of the turbulent spectra. As a 
result, significantly large differences in term of P(u>) 
are observed when using different turbulent spectra. 
This f eature is also well expl ained for Procyon (model 
C) in|Samadi et al. (2000b[). 



Besides a temperature effect, there is also a luminos- 
ity (or mass) effect although much smaller. This is 
illustrated with model F which is much hotter than 
model E but shows only a very weak sensitivity to 
the choice of the turbulent spectrum when comput- 
ing P(oj) (not shown here). In that sense, model F 
is closer to the cooler model E than to the other 
stars with similar effective temperatures. It is an 
intermediate case as it is less massive than for in- 
stance model C . Therefore the excitation region for 
this star extends deeper down compared to that of 
Procyon (model C). 

Figure |^ presents a star to star comparison of the 
oscillation power computed assuming the NKS. The 
plot is performed in the HR diagram in order to de- 
pict the dependence of the power spectrum P(lo) 
with the effective temperature and the luminosity. 
As the star becomes hotter and more luminous the 
frequency domain where P(oj) takes significant val- 
ues becomes smaller. This is due to the decrease of 
the cut-off frequency with increasing values of T e g 
and L. 



3.3. Observational constraints 



Observations of solar-like oscillations in Procyon 
have y et be en only clearly disco vered by Martic et al 



(1999 ) and Barban et al. (1999 ) in spectroscopi c sur- 



face velocity meas urements. However currently But i 
ler fc et al. (2000| ) claimed the detection of solar like 
oscillations in p Hydri. 

Most of the current ground based observations are 
mainly limited by the daily aliases. In particular 
they cannot yet provide the growth rates n which 
are necessary in order to compute P(oj) from Eq. ([!]). 
Thus we need observations of solar-like oscillations 
performed in such a way as to avoid the daily aliases. 
In particular space based experiments are particu- 
larly adapted for asteroseismology. The forthcoming 
space project COROT (Baglin & the Corot team, 
1998) based on photometr ic measurements will reach 
a noise leve l of 0.7 ppm (Auvergnc fc the COROT 



Team, 2000 ) and will thus detect oscillation ampli- 
tudes comparable to the solar ones (~ 2 ppm). Fur- 
thermore the instrument will continuously monitor 
several stars during ~ 150 days giving a frequency 
precision of ~ 0.1 /iHz such that accurate measure- 
ments of n will be available. 

Will this high accuracy enable us to constrain the 
theory of the stochastic excitation and the stellar 
turbulence ? The answer depends on the detection 
threshold and the accuracy in term of P(oj). For our 
stars here, n are not currently available. It is possible 
however to evaluate crudely a detection threshold for 
P(uj). Indeed, the power is related to the root mean- 
square of the surface velocity v s according to Eq. (|J) 
where r s is set to the radius at which T = T c ff. In 
the adiabatic assumption v s is in turn simply related 
to the luminosity fluctuation SL approximatively as 
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Figure 3. Oscillation power P versus frequency (v = lo/2tt) assuming the different spectra of Fig.l for the 
stars : A, B, C and E. Vertical error bars AP (plotted on the solid line) were obtained from Eq. mj assuming 
an accuracy for 77/ '2ir (resp. SL/L) o/O.l/xPz (resp. 0.7 ppm). For each star the detection threshold calculated 
according to Eq.(7) is represented by a dotted line. 



(Kjeldsen fc Bedding, 199E) 



SL/L oc v s T eff 1/2 



(4) 



Eq. (|^) and Eq. (|j) enable us to derive an approxima- 
tive relation 



P 



5L/L 



P Q \6L q /L q J (T cSV ICr 2 ) 



(5) 



where quantities with a are relative to the Sun. 
poldrcich fc Kumar (1991 ) derived a simplified equa- 
tion for 77 wh ich allows us to express 77/77© as (see 
Houdek, 1996|, section 3.4.4) 



77 _ L (H cj 2 lo Rf I- 1 
^to " Lq (Hcs^R^Iq 1 



(6) 



where c s and H are the sound speed and the pressure 
scale height respectively evaluated at the top of the 
convection zone and R is the star radius. Let s be 
the detection threshold in terms of relative luminos- 
ity fluctuation 5L/L expressed in unity of (5L/L)q. 
The detection threshold P ih (u)) in terms of oscilla- 
tion power can be expressed with the help of Eqs( 



as 



P th {uj) = b s 2 P™* 
with b = T cS L H 2 c- 4 I (T cff L H 2 c~ A 



R 2 C 2 M u; 2 

Rq £r,0 ( w max) w : 



2 

max 



(7) 



S JQ 



where w max and P™ ax are the frequency position and 
the power at the maximum value of SLq/Lq. In 
addition from Eq.(^|) relative error bar in term of P 
can be expressed as 



AP A5L Ari 

= 2— 1 - 

P 5L 7] 



(8) 



where 77 and SL are evaluated according to Eqs.(||||). 

In Fig.|| the detection threshold P th (uo) is plotted 
for the case s = 1 which corresponds to 3 times the 
noise level (rms) of the COROT instrument. In ad- 
dition the error bars AP according to Eq.(||) is plot- 
ted assuming COROT's performances, i.e. Av/2-k ~ 
0.1 /jHz and A5L ~ 0.7 ppm. Note that, as for the 
Sun, AP is larger at high and low frequency because 
SL becomes very smaller. 

We conclude that the differences between the spectra 
for the hotter stars (models A and C) are sufficiently 



■5 




Figure 4- Positions of the selected stars in the HR diagram (filled squared). For each star the computed power 
spectrum assuming the NKS is shown in an associated box. Abscissa extend from to the acoustic cut off 
frequency : the frequency range of the oscillation power narrows for more massive stars. 



large compared to error bars to identify the best tur- 
bulent spectrum. 



CONCLUSION 

This work shows that for stars hotter and more mas- 
sive than the Sun, the formulation of the power P{ui) 
injected into the solar-like oscillations is very sensi- 
tive to the way the turbulent stellar spectrum is mod- 
eled. For the Sun it is presently possible to constrain 
the solar turbulent spectrum from observations of the 
granulation. Such information is not directly avail- 
able for other stars. Evaluation of P(lu) derived from 
the seismic observations of the solar-like oscillations 
(amplitudes and damping rates) will therefore pro- 
vide constraints on the stellar turbulence. 

In addition for a given spectrum it is found that the 
shape of P(co) versus uj can change quite significantly 
from one star to the other. These differences are re- 
lated to the size of the convection zone which de- 
pends on the fundamental parameters of the stars. 
A star to star comparison of the power derived from 
observations will thus provide additional constrains. 



The accuracy of an observational determination of 
P(u)) depends on the accuracy of the measurements 
of the associated oscillation fluctuations 5L and the 
oscillations damping rates r\. A crude relation de- 
rived in the present work allows to evaluate the 
COROT detection threshold and the accuracy in 
terms of P(uj). It suggests that futur space observa- 
tions of such solar-like oscillating stars will measure 
acoustic powers P(uj) with an accuracy which will 
able us to identify the turbulent spectrum closest to 
reality. 
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